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We search for acoustic analogues of a spherical symmetric black hole with a pointlike source. We 
show that the gravitational system has a dynamical counterpart in the constrained, steady motion 
of a fluid with a planar source. The equations governing the dynamics of the gravitational system 
can be exactly mapped in those governing the motion of the fluid. The different meaning that 
singularities and sources have in fluid dynamics and in general relativity is also discussed. Whereas 
in the latter a pointlike source is always associated with a (curvature) singularity in the former the 
presence of sources does not necessarily imply divergences of the fields. 

<n : 
a.: 

<^ ■ I. INTRODUCTION 

One of the most intriguing features of classical general relativity is the presence of curvature singularities. The 
classical gravitational dynamics described by general relativity generates singularities both in the collapse of massive 
bodies and in the cosmological evolution of the universe. The physical meaning of these spacetime singularities has 
been widely discussed since their discovery. By now it has become conventional wisdom the point of view that the 
resolution of the singularities of general relativity has to be found in a quantum theory of gravity governing its short 
distance behavior. For instance, string theory allows for the quantum resolution of some singularities, but we are 
. far away from having a detailed description of how the Schwarzschild black hole (or the cosmological) singularity is 
£^ ' resolved. 

Building on a proposal of Unruh Q] , in recent years there has been a flurry of activity about the use of condensed 
matter systems to mimic various kinematical aspects of general relativity [2, 0, IH [a, 01 ■ Condensed matter systems 
provide us with analogue models of gravity that not only reproduce typical gravitational structures (event and cosmo- 
logical horizons) but can be also used to indicate an answer to fundamental questions of gravitational physics (field 
5^ | theory in curved spacetime and Hawking radiation) 0, la, 0. HH fl2l fl3L [Tj, [lg, fl7l Il8j . The hope behind 
bJ[). these efforts is that in a near future we could experimentally test the properties of the gravitational interaction using 
condensed matter systems analogues. 

One is therefore tempted to use analogue systems of gravity, in particular fluids, to try to mimic and understand 
spacetime singularities. Till a few months ago this idea would have been immediately rejected because the analogy 
gravity/fluids was formulated only at a kinematical level. Spacetime singularities in general relativity have a dynamical 
origin. In order to describe them using analogue models one needs to formulate the correspondence at a dynamical 
level. 

For spherically symmetric black holes such a dynamical correspondence has been recently found [19|. After gauge 
fixing, the four-dimensional dynamics governing a spherically symmetric black hole was shown to be equivalent to 
that governing the steady, constrained, cylindrically symmetric motion of a fluid. The results of Ref. fl9| provide us 
with a framework to describe black hole singularities using acoustic analogues. 

In this paper, extending the results of Ref. [l!j we will find the dynamical acoustic analogue of a spherical symmetric 
black hole with a pointlike source. We will show that the gravitational system has a dynamical counterpart in the 
constrained, steady motion of a fluid with a source (Sect. II). Although the discussion will be purely classical it will 
shed light on the different meaning that singularities and sources have in fluid dynamics and in general relativity. 
Whereas in the latter a pointlike source is always associated with a (curvature) singularity, in the former the presence 
of sources does not necessarily imply divergences of the fields. 

We will also discuss the implications of our results for black hole thermodynamics (Sect. Ill) and solve the constraint 
for the fluid dynamics in the case of a fluid with a generic power-law equation of state (Sect. IV). 
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II. GRAVITATIONAL DYNAMICS, POINTLIKE SOURCES AND FLUID DYNAMICS 

We start from generic four-dimensional Einstein gravity coupled to matter fields and a pointlike source of mass m. 
The action is (we adopt the notations of Ref. |20| and use natural units c = h = fcs = 1) 



-.</ j < l>j '\ -.'/ £.\u • - " ! j "i \/-g»j ^ , (i) 

where G is the Newton constant and Cmf is the Lagrangian for the matter fields. Cmf not only describes matter 
fields such as gauge and scalar fields, but may also contain a cosmological constant A. We shall consider for simplicity 
only the case in which the classical solutions for the matter fields are (or can be treated as) constant background 
configurations. In general the field equations stemming from the action will admit spherically symmetric black 
hole solutions, whose mass M will be equal to the mass m of the pointlike source. Moreover, these solutions will 
exhibit curvature singularities at the location of the source. The simplest case is represented by the Schwarzschild 
solution, obtained setting Cmf — in Eq. flip. Introducing a cosmological constant term or an electromagnetic field 
in we have the Schwarzschild-anti de Sitter or the Reissner-Nordstrom solution, respectively. 

We want to find a classical, non-relativistic fluid whose dynamics is equivalent to that governing the static, spheri- 
cally symmetric solutions of the theory When the pointlike source is not present, this problem can be drastically 
simplified by noticing that once spherical symmetry is imposed the gravitational theory Q becomes essentially two- 
dimensional (2D). Recently this result has been extended also to the case where sources are present [2l|. The 2D 
model is obtained by retaining only the radial modes in the action Q and it has the form of a 2D dilaton gravity model 
(see for instance Ref. |22|). A scalar field dj (the dilaton) parametrizes the volume of the transverse 2-dimensional 
sphere, 

ds 2 (i) =ds\ 2) + ^(j>dnl, (2) 



where A is a parameter with dimensions of length , related to the 4D Newton constant, G = A 2 . Inserting the 
ansatz © into the action Q and performing the Weyl rescaling of the 2D metric 

9, (3) 



needed to get rid of the kinetic terms for the dilaton, one obtains a 2D model characterized by a dilaton potential 
V(cb) and a coupling function W((f>) 21], 



A = A G +A M = \J d 2 x^—g{4>R + \ 2 V{i,))-m J dtW^-g^^-. (4) 

The dilaton potential V(4>) depends on the form of the matter Lagrangian Cmf, whereas the coupling function W{4>) 
reads 

W{<t>) = (2d J )- 1/i . (5) 

For static solutions in the Schwarzschild gauge, 

ds 2 = -U(r)dr 2 + U^irjdr 2 , (6) 
with a source particle at rest in the origin, the field equations coming from the 2D action (@J become plj . 

d 2 U ^dV _ dW 



dr 2 



= A — — - 2m— S(r), 



dr dr 

d 2 <b dU deb , 
2U—^ + - — i~ = X V - 2mW6(r). 
dr z dr dr 

In order to find the acoustic analog of the gravitational dynamics Q we will consider separately the gravitational 
dynamics in vacuo, m — 0, and in the presence of the source m/0. 
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A. Gravitational dynamics in vacuo and fluid dynamics 

In this case the problem of finding a fluid dynamics equivalent to the gravitational dynamics J7J) has been already 
solved in Ref. |l9j |. However, in that paper the dilaton has not been considered as dynamical, but rather as constrained 
by the equations of motion to be proportional to the spacelike coordinate r of the 2D spacetime. Here we will improve 
the derivation of Ref. U| by considering the full gravitational dynamics with both degrees of freedom U and 0. 

When considering the field equations J7J) in the vacuum or, equivalently, the spacetime region r ^ away from the 
source, we can define a scalar mass function M(r) j23| 

M(r) = y( X2 J V ^ d<t> ~ ( V( ^ 2 ) ' ^ 

which is constant by virtue of the field equations Q, with Fq arbitrary normalization constant. On shell we therefore 
have M(r) = M = const, where M is the mass associated with the classical solution of the field equations under 
consideration. 

Using the mass function M(r) we can find an equivalent form for the system Q. Only two of the three equations 
in J7J are independent and the the system can be written, for <f> ^ const, as 

f^=AV, ^ = 0. (9) 
dr dr dr 

By subtracting the third from the second equation we can also rewrite Eqs. JJJ in a second equivalent form, 



^ = A 2 V, ^ = 0. (10) 
dr dr dr 1 

The existence of two equivalent forms for the equations of the gravitational dynamics can be traced back to the 
existence of two constant of motion: M and (j>o = (l/\)d<fi/dr. In the solution of Eqs. Q) these constants of motion 
appear as integration constants, 

-2 1 j, , 2M( t ) o\ , 2 , t(X\ 2M <Po\^ 1 , 2 



= -tf UW) - dr' + ^ IJM - 1 dr 2 , = o Ar, (11) 

where J = J Vd<j>. Moreover, they are related one to the other in a geometrical way: d(p/dr determines the Killing 
vector of the metric <|11[1 = ^"dvfy, whereas M is its associated conserved charge. In general the solution <|llf) 
describes a black hole with an event horizon at r = r^, with J(r^) = 2Mcj>o/\. 

The form @ of the gravitational field equations appears more suitable for making contact with fluid dynamics. For 
this reason in the following we will consider the equations for the gravitational dynamics in the form given by Eqs. 
@. With this choice the constant <j>o becomes irrelevant. Therefore we will set <fio = 1, and correspondingly Fq = 1/A 
in Eq. ©. 

The gravitational metric JTJ can be related to the acoustic metric associated with the steady, cylindrically sym- 
metric flow of a barotropic, inviscid and locally irrotational 3D fluid [l9| . 



ds 2 = ^ [- (c 2 - Vq) dt 2 - 2v dxdt + dx 2 ] , (12) 



where x is the coordinate along the flux tube, po(x) is the density of the fluid, vq(x) its velocity and c[x) — \J dP/dp 
is the velocity of sound (P is the pressure of the fluid). Owing to the cylindrical symmetry of the flux tube the 
fluid parameters p, vq , c, the section of the flux tube and the potential for the external forces acting on the fluid 
are independent of the transverse coordinates y, z of the flux tube. For notational convenience, in l|12|) we use the 
dimensionless variable po = po/A 4 . The mapping between the gravitational (|11|) and acoustic 112|) metric is realized 
by the coordinate transformation 

Podx, r = t+ / dx^ — — o, 



2M _ po 



[c 2 -v 2 ). (13) 
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Considering a flux tube of section A(x), the dynamical equations governing the motion of the fluid are the Euler and 
the continuity equations 

dvo dP dip 
Pqv — + — + p Q — = 0, 

dx dx dx 

dT d , T , , 
-=-(- Po v A) = 0, (14) 

where T is the flux of fluid mass and ip is the potential for external forces acting on the fluid. 

Let us now compare the gravitational equations with the fluid-dynamical equations obtained defining the variables 

X = ^(c 2 -vZ), Y = Po c, F = ln(^-). (15) 
c \PoJ 

The Euler and continuity equations take, respectively, the form 

dX _ ^dY x dF 2e _F# 
dr dr dr dr ' 



— y/A*Y(Y-X) = 0, (16) 
dr 

where we have introduced the dimensionless flux tube section A = X 2 A. 

The fluid dynamic equations 1(161) are put in the gravitational form @ by the identification 

X = U, y/A^Y(Y-X) = ^-, (17) 

(actually, on the right hand side of the second equation one could choose any function of M(r)), and by introducing 
the constraint 

„ dY dF p dib , 9 . , dr , . 

2 ^~ X ^ + 2e Tr = XV ^Tf (18) 

In Eq. 1(17(1 the A factor has been introduced for dimensional reasons. Notice that the identification (|17|l and the 
constraint (|18|) hold both for a flux tube of constant and non-constant section. 

In terms of the solutions of the gravitational field equations Ijllf) the identifications (|17|l read 

X=(,J(\r)- 2 -fy ^ A *Y(Y-X)^^ = ^, (19) 

where now M is M(r) evaluated on-shcll. The identification given by the second equation in (|19f) is very natural and 
has a simple physical meaning. It identifies the conserved quantity of fluid dynamics, the flux of matter fluid !F, with 
the conserved quantity of the gravitational dynamics, the black hole mass M. 

B. Gravitational dynamics with pointlike sources and fluid dynamics 

Let us now generalize the previous results to the case m/0. Also in presence of a pointlike source the equations 
can be rewritten in terms of M (r) , in a form similar to , as 

M4£ =x * Vi d M™ W d A 6{r) . (20 ) 
dr dr dr A dr 

Introducing the variables (|15JI we find that also in this case the identification 1)17(1 . allows us to find the acoustic 
analog of the gravitational system 120(1 . The Euler equation in 1(161) and the constraint ((18(1 are unchanged by the 
coupling to matter, and are still equivalent to the first equation in ((20(1 . The second equation in 1(20(1 implies a 
modification of the continuity equation in l(16|) . which acquires a term proportional to the parameter m, 

UA P ov) = ^W d ^(r). (21) 
dr \ z dr 
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This can be recognized as the continuity equation in the presence of a flux source of strength 

rn 



$ = 
A 2 



S(r). (22) 



=o 



The acoustic analog of the gravitational singularity is therefore simply a source term for the flux of the fluid. Owing 
to the cylindrical symmetry of our flux tube the source is planar. Using the first equation in ((13(1 one can trade in 
Eqs. i(21|) . 1(22(1 the radial coordinate r for the the flux tube coordinate x. In terms of the coordinate x the source is 
located at the position x = xq corresponding to r — and it extends along the perpendicular y, z directions to cover 
the whole section of the flux tube. 

In the gravitational description the position of the source is associated with a singularity of the dynamics (curvature 
singularity). This is not the case with the acoustic description where the delta function singularity at r — appears 
as a pure source term for the fluid flux. The Euler equation and the constraint 1(18(1 are unaffected by the presence of 
the source. We see here an important difference between the gravitational and the acoustic case. In general relativity 
the presence of a pointlike source implies a (curvature) singularity for the gravitational field. This is not the case 
for the fluid where the presence of the source is completely disentangled from the appearance of a singularity in the 
dynamics. 

The solutions of the gravitational field equations J7J have been derived in Ref . j2lj . The delta function singularity 
at r — is generated by taking the solution as function of \r\, 

tf = M_ 7j 4> = a\\r\+(3, (23) 
cr- 

where as usual J = J d<pV and a, (3, 7 are integration constants satisfying 

A third equation constraining the values of the integration constants has to be added to ensure the equality between 
the mass M of the solution and the mass m of the source |2l|): 

2m 

^=A^ (25) 

Notice that, as expected, solution ((23(1 differs from the vacuum solution (jl 1|) only at r = 0. For r > the two solutions 
are related by a translation of the coordinate r and a scale transformation of r and r. One can now easily work out 
the coordinate transformations relating the gravitational black hole (|23|l to the acoustic black hole 112(1. They are 
given by Eq. 1(13(1 with the third equation replaced by 

a- 2 J- 7 =^(c 2 -, 2 ). (26) 



III. BLACK HOLE THERMODYNAMICS 



The correspondence between gravitational and fluid dynamics allows us to straightforwardly define the thermody- 
namical parameters, T a , M a , S a associated with the acoustic black hole, which satisfy the first principle dM a = T a dS a 
[T^ | . The relevant formulas expressing the thermodynamical parameters in terms of the variables X, Y, F of Eq. 1(15(1 
have been already given in Ref. 19] . However, the acoustic black hole thermodynamics developed in that paper is 
not completely satisfactory. The thermodynamical parameters associated with the acoustic black hole are a simple 
"translation" of those associated with the gravitational black hole. A physical interpretation of them, in particular 
of the mass M a = M, in terms of fluid parameters is difficult. 

The results derived in the previous sections of this paper, in particular Eq. I(19() . allow us to give a more transparent 
"fluid-dynamical" interpretation of the thermodynamical parameters associated with the acoustic black hole. From 
Eq. 1(19(1 it follows that the black hole mass M is the flux of matter fluid measured in units of A, M a = T j\. Because 
T (and M) are constant of motion we can evaluate them on the horizon r = r^. We have, 

M a = M(r h ) = XA(r h )Y(r h) . (27) 

The entropy of the acoustic black hole is — 27tAtv For a flux tube of constant section the entropy is 

proportional to the total mass of the fluid M.{rh) inside the horizon. In fact from Eq. I(13|) it follows 

f rh 2?r 

S a = 2n\ J Po dx = jjM(r h ). (28) 
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The thermodynamical parameters satisfy the first principle. Making use of Eq. 1(27(1 we therefore have for the black 
hole temperature, 

dM a 1 d 1 dT(r h ) 

Ta = ^ = 2^ M^rs)] = 2^^^- (29) 

The temperature of the acoustic black hole measures the rate of change of the flux of fluid mass when the position 
of the horizon is changed. For a flux tube of constant section, for which is proportional to the mass of the fluid 
inside the black hole, the temperature measures the rate of change of the flux when the mass of the fluid inside the 
hole is changed. 

IV. SOLUTIONS OF THE CONSTRAINED FLUID DYNAMICS 

In the previous sections we have shown that the gravitational dynamics governing the spherically symmetric solution 
of the action Q has an acoustic analogue given by a fluid dynamics constrained by Eq. 1(18(1 . If the geometry of the 
flux tube and the external forces are fixed, this constraint determines the equation of state of the fluid. The physically 
relevant situation, which we consider here, is the opposite. One works with a fluid with a given equation of state and 
Eq. HI 8(1 becomes a constraint on the geometry of the flux tube and/or on the external forces acting on the fluid. 

In Ref. the constraint ((18(1 has been solved for a perfect fluid. In the following we will solve the constrained 
fluid dynamics for a fluid with a generic power-law equation of state 

P = ^-Po- (30) 
n 

where a and n ^ are real dimensionless constants. This equation of state describes almost all physically interesting 
fluids: perfect fluid (n = 1), Bose-Einstein condensate (n — 2), Chaplygin gas (n = — 1), fluid mechanics in d spatial 
dimensions invariant under the nonrelativistic conformal group SO(l, 2) (n = 1 + 2/d) [24j. 
In terms of the variables ((15(1 the equation of state reads 

y3-« e (n+l)F = a ^ (31) 

The particular cases n — 1 (perfect fluid), n = — 1 (Chaplygin gas) n = 3 ( conformal invariant fluid mechanics in 
d = 1 spatial dimensions), corresponding respectively to Ye F — const, Y — const, F = const will be considered 
separately. 

Using Eq. 1(31(1 in (the inverse of) Eqs. 1(15(1 we find the fluid velocity, the fluid density and the speed of sound 



u = a— ywi) y/Y-X, (32) 

p = a-^TLY^TL, c = a^Y^i. (33) 
We will discuss separately the two cases of a flux tube of constant and non-constant section. 

A. Flux tube with constant section 

In this case A = const and we can use the continuity equation to solve for X — X{Y): 

a 2 , s 

X = Y- Y , (34) 

where a = M /(AX). Using Eq. ((34|) into Eq. IMl'l) we can express also vq as a function of Y: 

v o = a (t?) " +1 ( 35 ) 



The constraint 1(18(1 can be written as 



Y 



dX n , v dF n F dip 



dY 2 + X dY- 2e d^ = °> ^ 
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which using Eq. (|31|l can be easily solved for ip(Y), 

$ = - a ^T ( ^y-^r + -J-Y 2 ^) . (37) 
\ 2 n — 1 / 

The acoustic horizon forms at Y — a, corresponding to X = 0. The subsonic region (X > 0), describing the external 
region of the gravitational black hole, is given by Y > a. The source, corresponding to the black hole singularity, is 
located at Y = Y s with < Y s < a. The fluid parameters vo,po,c and the external potential ip remain finite both 
on the horizon and on the location of the source. In the case under consideration the constraint necessary to have a 
correspondence between gravitational and fluid dynamics takes the form of a constraint on the form of the external 
forces acting on the fluid, given by Eq. I|37[l . Moreover, this constraint implies a "null force condition" on the acoustic 
horizon. In fact the external potential ip given by Eq. I|37|l becomes extremal for Y = a (a local maximum for n > — 1 
and a local minimum for for n < —1). 

Once a specific action for the four-dimensional gravity model QJ, and hence a specific V(tfi) in the 2D model 
is given, one can use Eq. (|34|l to find Y(r). Use of Eqs. (|33|l . Ij35|l and l|37(l allows one to get the fluid parameters as 
a function of r. 



B. Flux tube with non-constant section 



In this case A is not constant and we can use the continuity equation to express it as function of X and Y, 

* = Yt&ry < 38 > 

The presence of external forces is now an unnecessary complication, and therefore we set them to zero. With ip = const, 
and using Eq. (|3*T)l the constraint (|T%|) becomes 

dX n+lYdX~ ' (dyj 

which is readily solved to give 



l(Y-u 2 Y-^, (40) 



X = 

where cu is an integration constant. Inserting the previous equation into Eqs. (|32|l . I|33|l and l|38|) . we get the solution 
of the constrained fluid dynamics. The acoustic horizon forms at Y = Yh = and the source is located at 

Y = Y s , with < Y s < Yh. The supersonic region, corresponding to the black hole interior, is given by Y s < Y < Yh- 
Conversely, the subsonic region (black hole exterior) is given by Y > Yh- The fluid parameters Vq , po and the speed 
of sound c remain finite both at the horizon and at the location of the source. 

We can also show that for n > — 1 the acoustic horizon always forms at a minimum of the section A(x), i.e. the 
flux tube must have the form of a so-called Laval nozzle. In fact it follows from Eq. Q38J1 and Eq. Q39JI that 



dA AX 1 



dX Y - X 2(n + 1)Y - (n - 3)X 



(41) 



Taking into account that by definition Y > and that Eq. (|38|l requires Y > X, one finds that for n > — 1, 
dA/dX = at the horizon (X = 0), whereas it is positive (negative) for X > (X < 0). This is a highly non- 
trivial result. The constraint, which is necessary to have a correspondence between gravitational and fluid dynamics 
becomes a geometrical constraint on the form of the flux tube. This geometrical constraint forces the fluid to develop 
an acoustic horizon. 



C. Particular cases: n = 1, 1,3 

The case n = 1 corresponding to a perfect fluid has been already discussed in Ref. 0. For n — — 1 we have the 
Chaplygin gas. The pressure is negative and the equation of state becomes P — —\ 4 a 2 /p , which in terms of the 
new variables l|15f) reads Y = a. In the case of a flux tube of constant section, constant Y implies, owing to the 
continuity equation, also constant X. This forbids the realization of an acoustic Chaplygin analogue of a gravitational 
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black hole, independently of the form of the potential V. This difficulty can be circumvented by considering a flux 
tube of non-constant section. In this case, considering for simplicity ip = 0, the constraint (|18|l can be solved to give 
exp(— F) = bX, with integration constant b > 0. The fluid parameters take therefore the form 



This solution is defined only for X < a and becomes singular on the black hole horizon X = 0. 

For n = 3 the equation of state implies F — const and the solutions take a very simple form. In the case of a flux 
tube of constant section Eqs. lO > E3J 7 CH3>! ijSZjl gi ye 

-o=a(^) 1/2 , p =(^f 1/2 , c^^F) 1 / 2 , (43) 

y = i(x-v/^T4^), ^ = -|^ + y^. (44) 

In the case of a non-constant flux tube section we get instead 

v = \ja ( d- ^- ), po = \/^ ( d+ 77) c = ap , A=-^-. (45) 

where <i > is an integration constant related to the integration constant lu appearing in Eq. I|4U|I . It is evident from 
the previous formulae that —2c? < X < 2d. The horizon is located at X — 0, which is also a minimum of A. The 
section of the flux tube diverges at X = ±2d. Therefore the acoustic black hole cannot describe the entire black hole 
spacetime. 





V. CONCLUSIONS 



The resolution of the singularities of general relativity relies heavily on the understanding of the short distance, 
quantum behavior of gravity. Nevertheless, investigations of singularities at the classical level may also help to improve 
our understanding of the subject. In this paper we have presented acoustic analogues of the black hole singularities 
generated by a pointlike source. The acoustic analogues of the singularities are constructed finding a fluid whose 
dynamics is completely equivalent to that governing the spherically symmetric solutions of (gauge fixed) general 
relativity with pointlike sources. 

This result enhances the analogy between gravity and condensed matter systems. This analogy is not restricted to 
kinematical features of gravity such as acceleration horizons. We can use a fluid to mimic the full black hole dynamics, 
spacetime singularity included. Analogue models of gravity have also the nice feature of being, at least in principle, 
experimentally testable in laboratory. In the near future this could open the way to experimental tests of the full 
classical black hole structure realized using condensed matter systems. 

Another important result of our investigation concerns the notion of singularity and source in classical field theory 
and, in particular, in general relativity In a classical field theory, such as Maxwell electromagnetism, the (delta 
function) singularity associated with a pointlike source is conceptually independent from the field singularity, the 
point where the fields diverges. It is the classical dynamics (e.g. the Laplace equation) that forces the identification 
of the two types of singularities. The delta function singularity becomes the source of the Green function for the 
field. In general relativity this relationship between pointlike sources and field singularities becomes even deeper. The 
equivalence principle treats on equal footing sources and test particles, and the field singularities become geometrical 
(curvature singularities). 

In this paper we have shown that for a classical field theory describing the motion of a fluid we may have a slightly 
different situation. The presence of a source for the mass of the fluid is not necessarily related to a divergence of the 
fields describing the fluid. It is related with a much milder form of singularity, a cusp point in the field configuration. 
This result is consequence of two features of the the fluid dynamics we have considered. First the effective theory which 
describes the fluid is two-dimensional. In two dimensions the delta-like singularity of the source can be generated by 
a second derivative of |r| in the Green equation. Second, differently from other field theories, in fluid dynamics the 
source couples, through the continuity equation, to the flux of matter fluid only. 

By mapping gravitational and fluid dynamics, we have also shown that the singularity of spherically symmetric 
black hole spacetime generated by a pointlike source can be transformed in the milder singularity of fluid dynamics 
described above. The responsible for this miracle is the 2D dilaton-dependent Weyl rescaling of the metric @. The 
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price that we have to pay for this is that the Weyl transformation itself is singular at the position of the source (see 
eq. (J23J). The Weyl transformation changes the the 4D black hole singularity into a 2D cusp singularity. 

In this paper we have considered only spherically symmetric solutions of general relativity. It is an interesting open 
question if our results hold also in more general situations. It seems very difficult to mimic a gravitational dynamics 
with propagating degrees of freedom or a system with two or more interacting masses, using a fluid. Nevertheless a 
description similar to that presented in this paper may still hold for classes of solutions with high degree of symmetry, 
such as for instance cosmological solutions. 
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